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The interpretation of y(4140) as a D^Dl molecule is studied dynamically in the one boson 
exchange approach, where cr, r) and exchange are included. Ten allowed D*D* states with low spin 
parity are considered, we find that the J^*^ = 0"''^, 1^ , 0~^, 2^"'" and 1~~ D^Dl configurations 
are most tightly bound. We suggest the most favorable quantum numbers are J^"^ = 0"*"+ for 
y(4140) as a D*D* molecule, however, J^'~^ = and 2"*"+ can not be excluded. We propose to 
search for the and 1 partners in the J/tpr] and J/ijjri' final states, which is an important test 
• of the molecular hypothesis of y(4140) and the reasonability of our model. The 0''"'' B*B* molecule 

' is deeply bound, experimental search in the T{lS)(j> channel at Tevatron and LHC is suggested. 

o : 

CN ' PACS numbers: 

D ■ 

c/^ : 

CN ' I. INTRODUCTION 

^ . 

I— I. Recently the CDF Collaboration reported a narrow structure F(4140) in the decay 5+ X+F(4140) followed by 
i-C ' y(4140) — > J/ipcj) with a statistical significance of 3.8 a. The mass and the width of this structure are fitted to be 
P^'. 4143.0 ± 2.9(stat) ± 1.2(syst) MeV and 11.7i5 o(stat) ± 3.7(syst) MeV respectively [10. Since the quantum numbers 

of both J/tp and (j) are l'^{J^'~^) = 0~(1~~), F(4140) is an isospin singlet with positive C parity. 
D , F(4140) is very similar to the charmonium like state y(3940), which is observed in B ^ KY{2>%AQ) KJ/ipu 
, ■ [1, It was argued by the CDF Collaboration that y(4140) can not be a conventional charmonium state, because 
\ a cc charmonium state with mass about 4143 MeV would dominantly decay into open charm pairs, and the branch 
PsJ ■ ratio into the double OZI forbidden modes J/ip4> or J/tpui is negligible. Although the decay F(4140) J/i^(t> 
K*" ' is unusual, the decay Xbi 2(2P) ijjT{lS) has been observed Therefore F(4140) as a conventional cc state 
can not be ruled out completely due to the scarcity of current experiment data. Comparing with the theoretical 
predictions for the charmonium spectrum [1, 0|, we suggest that F(4140) would most likely to be the 2^D2 state 
with I^{ J^^) = 0+(2~+), if it is a cc charmonium state. It would dominantly decay into DD* , D*D* and DgD* in 
• ■ this scenario. A possible explanation for the unusual decay mode y(4140) J/ip(f> is that the rescattering through 
^ ; DsD*jDlDs or may be_ responsible, i.e. F(4140) DsD*{D*Ds or D*D*) J/i^(j). Mixing between the 

■ charmonium state and the D*D* molecule may also contribute to this extraordinary decay. 

I There are already some theoretical discussions about this structure. The authors in Ref. Q suggested that F(4140) 
• • . is a molecular partner of y(3940), and that it is a J^'~^ = 0++ or 2++ D*D* molecule. In Ref. Q, the author argued 

' that F(4140) is a D*D* molecule or a exotic 1~+ charmonium hybrid. The interpretation of F(4140) as a 0++ Dt^t 
^ , molecule has been studied by QCD sum rules [l3, [HI, however, different conclusions were reached. The decay of 
^ ■ y(4140) as a hadronic molecule or conventional charmonium x'co i were discussed as well [l^ [isj. 

\ In the past years, some new states near threshold have been observed, and they are usually suggested to be hadronic 
molecules [H, [H, M, E.lil, [H, El. Being different from familiar hadr onic molecule candidates, F(4140) is about 
81.6 MeV below the D^D* threshold, its binding energy is not small if it is identified as a DID^ molecule. There 
are various theoretical methods to study the hadronic molecule dynamically, such as the QCD sum rule, the QCD 
effective field theory [2^ [2l|, the potential models with pairwise interactions between quarks fl6l . [TqI and so on. 
Tornqvist's original work on hadronic molecule from one pion exchange is especially attractive [22|, |23|. Guided by 
the binding of deuteron, he performed a systematic study of possible deuteronlike two mesons bound states with long 
distance one pion exchange. Being different from the above mentioned theoretical methods, he took into account 
the contribution of the tensor force. Since the tensor force turns out to be very important in the deuteron binding, 
one expects that it would make a significant contribution to the binding of the general hadronic molecule. At short 
distance, the interaction should be induced by the interactions among the quarks in principle. However, a detailed 
and reliable modeling of the short range interaction is not an easy matter, and various phenomenological models have 
been proposed [2^. [25j. although one pion exchange is expected to be dominant for the hadronic molecule. Inspired 
by the nucleon-nucleon interactions, we further extended the one pion exchange model to represent the short distance 
contributions by the heavier mesons rj, tr, p and to exchanges in Ref. [2^. From these studies, we learn that one 
pion exchange really dominates the physics of hadronic molecule, and the tensor forces indeed plays a critical role in 
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FIG. 1: The interaction potentials Vc{r), Vs{r), Vrir) and Vls('") with A = 1400 MeV. The sohd hne represents the central 
potential Vc{r), the dotted, dashed and dash-dotted lines denote Vs{r), Vrir) and Vls('') respectively. 

providing the binding, which results in the S-D wave mixing. Since in general the molecular state is weakly bound, 
the separation between the two hadrons in the molecule should be large. Consequently the dominance of one pion 
exchange can be easily understood. A remarkable success of the model is its prediction for the existence of the famous 
X(3872) long ago [H, [l^l- Based on one pion exchange, Frank Close et al. recently suggested that the charmed 
mesons molecular state is a possible solution to the enigmatic states observed above 4 GeV in e+e~ annihilation 
(27t . Therefore we think that the one boson exchange model can give a reasonable and reliable description for the 
hadronic molecule. In this work, we shall study the possible molecular states of the D*D* system and discuss the 
most favorable quantum numbers of y(4140) as a D*D* molecule in the one boson exchange model. Ten allowed 
D*D* states with total angular momentum smaller than 3 will be discussed. 

The paper is organized as follows. In section II. the effective potentials of the ten D*D* molecular states are given. 
The numerical results and the possible D^Dl molecules are presented in section III. Section IV discusses the most 
favorable quantum numbers of y(4140) and the implications, if it is interpreted as a D*D* molecule. Finally we give 
our conclusions and some discussions in section V. 

II. THE EFFECTIVE POTENTIALS FOR THE D'^D* STATES WITH LOW SPIN PARITY 

In the one boson exchange model [23,[2^,[2^ , the effective potentials between two hadrons are obtained by summing 
the interactions between light quarks or antiquarks via one boson exchange, and one boson exchange only between 
up quarks or down quarks has been considered. We note that the present one boson exchange is different from the 
well-known Goldstone boson exchange constituent quark model [28| , where one boson exchange between quarks inside 
a hadron is used to describe the light baryon spectrum. Since strange quark is involved in the present problem, we 
should extend our consideration to the one boson exchange between strange quarks. Because I?* is an isospin singlet, 
the effective potential between D* and D* is induced by rj, a and (j) exchanges, whereas the tt, p and lo exchanges 
don't contribute. Taking into account the overall sign difference (—1)'^ between the quark-quark interactions and the 
quark-antiquark interactions, where G is the G-parity of the exchanged meson, we obtain the effective potential 

V{r) = V„{r) + V^r) - V^{r) = Vc{r) + Vs{r) a, ■ a, + YT(r)S,,{i) + VlsWL • S,, (1) 

where Vniy), Va{r) and V0(r) respectively denote the effective potentials from ry, g and exchange between two 
quarks, and their expressions are given explicitly in Eq.(9)-Eq.(ll) of Ref. [1^. Noting that the exchange distin- 
guishes the D*Dl from the exotic D*D* systems. The subscript i and j are the indexes of the strange quark and 
strange antiquark. Sij{r) = 3(cr.i • f)(<Tj • f ) — cr^ • a-j is the well-known tensor operator, = ^{(Ti + (Tj) is the total 
spin of light quarks, and L is the relative angular momentum operator. Comparing with the general results in Ref. 
[201, the isospin dependent interactions vanish here. We would like to note that the spin operator CTi or CTj only acts 
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on the strange quark or antiquark here. The potentials Vc{r), Vs{r), Vrir) and Fls(^) £^rc given as follows: 

/ \ diss Tj \ diss "^l U f\ \ ^isS TT I K ^ , diss + '^d'pSsUsS '^^ „ , . s 

Vc{r) = — —maHa{Kma,r) - ^— ^— |iJi(A, m^, r) - m^Ho{A,m^,r) + — — ^--^7Ji(A, m^, r) 

= -1^12;^^^ + 6;;if^^(^'"^*'") 

-j^ / N 9ijss rj- / A \ 1 {9<i}SS + .ftfissY ^(p TJ f A \ 

^H0 = ^Y2;^i?3(A,m„r) + ^H,{A,m,,r) 

M I \ diss „ (K ^ , ^dlss + ^9<t>SsUsS „ . 

VLsir) = -17^^2(A,m.,r) + —H,iA,m,,r) (2) 

where rus is the mass of the constituent strange quark, and we take rris = 0.55 GeV in this work, m^, rria and are 
the masses of 77, tr and (j) mesons respectively. 5aiss(-^ = V, o' and 0) and /^ss denote the effective coupling constants 
between the exchanged meson and the strange quarks. In the above formulae, we have introduced form factor at each 
interaction vertex, and the form factor in momentum space is taken as 



where A is the so-called regularization parameter, m and q are the mass and the four momentum of the exchanged 
boson respectively. The form factor suppresses the contribution of high momentum, i.e. small distance. The presence 
of such a form factor is dictated by the extended structure of the hadrons. The parameter A, which governs the 
range of suppression, can be directly related to the hadron size which is approximately proportional to 1/A. However, 
since the question of hadron size is still very much open, the value of A is poorly known phcnomcnologically, and it 
is dependent on the models and applications. In the nucleon-nucleon interactions, the A in the range of 0.8-1.5 GcV 
has been used to fit the data. The extended structure of hadrons also has the following obvious consequence: because 
the mass of the exchanged meson determines the range of the corresponding contribution to the D*D* interactions, 
one should restrict oneself to meson exchange with the exchanged meson mass below a certain value, typically on the 
order of the regularization parameter A. The ry, a and (p exchange are considered in the present work, therefore the 
value of A should be at least larger than the 4> meson mass. The functions i/o(A,m, ?'), Hi{A,m,r), i?2(A,m,r) and 
H^^A, m, r) involved in Eq.Q are given by 

HoiA,m,r) = ^ (e"- - e"^'') - ^^^^ e^^^ 
mr 2mA 



H,iA,m,r) = -^(e"'"'^ - e'^O + ^^^^^^^^ e"^' 

1 X 1 ... 1 X A 1 



H,iA,m,r) = (l + _) —-e"""- - (1 + -) -— -e' 



Ar 

mr' rn?r'^ Ar mm'^r^ 2iri? mr 

3 3,1 ,,3 3 , A2 e-^- A^ 



i?3(A,m,r) = l + _+ _e-"-- 1 + — + — ^ — ^_ 1 + Ar (4) 





States 


0++ 


S'o, Do 








'Po 


1++ 


'D, 


1-+ 


'Pi 


2+- 




1"" 




2++ 


^1)2, 'S2, ''£'2, ''G2 


2-+ 


'P2, ^F2 


2"- 


'P2, 'P2 



TABLE I: The allowed states of the D^Dl system with low spin parity, where we cut off the total angular momentum to J = 2. 
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The four components of the interaction potentials Vc{r), Vs(r), Vrir) and VLsi^) ^^'^ displayed in Fig. [TJ We see 
that the central potential Vc{r) is negative, whereas the remaining three components are positive. These potentials 
are the cancellation results of different contributions, especially for the tensor potential Vrir), whose amplitude is 
smaller than that of the other three potentials. For the DID* system, the spatial parity is determined by P = (— 1)^ 
and the C-parity is C = (— l)^"'"'^, where L is the relative angular momentum between D* and D*, and S is the total 
spin of the system. We cut off the total angular momentum to J = 2, the allowed states with low spin parity are listed 
in Tabic [H The matrix elements of operators cr^ • crj, Sij(i) and L • can be calculated straightforwardly with the 
help of angular momentum algebra, and the results arc given analytically in the appendix of Ref. 26] . Consequently 
we can write out the one boson exchange potentials for the allowed states in matrix form as follows 
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(6) 



^0-+ W = Vc{r) - Vs{r) - 2VT{r) - VLs{r) 



(7) 



Vi++(r) = Vc{r) + Vs{r) - Vrir) - -^lsW 



Vi-+(r) = Vc{r) - Vsir) + V^r) - -V^sir) 



V2+- (r) - Voir) - Vsir) + Vrir) - -14sW 



(8) 
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In the following, wc will perform the same analysis as for the deuteron and the possible heavy flavor molecules in 
Ref. [26| . One can then determine for which quantum numbers the one boson exchange potential is attractive and 
strong enough so that D*Dl bound states are expected. 
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III. NUMERICAL RESULTS AND POSSIBLE D*,Dl MOLECULES 

The input parameters in our model are the involved meson masses, the effective coupling constants gMss{M = 77, 
a and </>) and /^ss, and the regularization parameter A. The meson masses are taken from the compilation of the 
Particle Data Group [H: = 547.853 MeV, m„ = 600 McV, = 1019.455 MeV and mr,. = 2112.3 MeV. In the 
SU(3) flavor symmetry limit, the coupling constants gMss(J^ ~ ih and and f^ss between the exchanged mesons 
and strange quarks are related to coupling constants gMqqi-^ ~ and </>) and f^qq between the exchanged mesons 
and up/down quarks via the following relations 

9i)ss '^g^iqqi Scrss Saqqi 9(pss — V^.9cjqq; feftss — ^^/^ fujqq (^^) 

The coupling constants gMqqi-M ~ rj, a and uj) and fu}qq can be estimated from the phenomenologically known rjNN, 
aNN and ojNN coupling constants. Riska and Brown have explicitly demonstrated that the nucleon resonance 
transition couplings to tt, p and oj can be derived in the single quark operator approximation [30j . which are in good 
agreement with the experimental data. Adopting the same method, we can straightforwardly derive the following 
relations between the meson-quark couplings and the meson- nucleon couplings p6l.l30j. 

mq 1 1 r rriq 1 niq 

9mq — 9vNN, gaqq — -;;9crNN, 9ujqq = 7; 9i^NN, Jujqq — JujNN ~ Ki; )9uiNN i-Loj 

niN 6 S niN 3 TTiN 

where mjv is the nucleon mass. In this work, the effective couplings between the exchanged bosons and the nuclcons 
are taken from from the well-known Bonn model [3l|: 5^jv7v/(4''') = 3.0, g^NN/i^^) — 7.78, 5'5AfJv/(4''') = 20.0 and 
f^NN/i^^) ~ ^- nature and in QCD, the flavor SU(3) symmetry is broken by non-equal masses of the up and 
down quarks and the strange quark or the electromagnetic effects. It is commonly believed that the error of SU(3) 
predictions is approximately 20% — 30%. Consequently the uncertainty of the coupling constants 9jjssi 9a-ssi 9ct>ss and 
f^ss is about 20% — 30% as well. As a demonstration for the consequence induced by the uncertainties of the effective 
couplings, all the coupling constants would be reduced by 20 percents later, and the corresponding predictions are 
analyzed seriously. Taking into account the centrifugal barrier and solving the coupled channel Schrodinger equations 
numerically, then we can obtain the predictions for the binding energy and the static properties, which are listed in 
the tables of the Appendix. We notice that these predictions are rather sensitive to the regularization parameter A 
and the effective couplings, this is common to the one boson exchange model [1, [2^, [33|. We also find that the binding 
energy increases with A, the reason is that increasing A increases the strength of the potential at short distance. 

For the 0++ D*Dl state, the system can be in ^5*0 or ^Dq configuration, this is very similar to the deuteron, 
which can be in ^5*1 or ^Di. The S wave state mixes with the D wave state under the tensor force, as is shown 
explicitly in Eq. ([5]) . The energy of the system would be lowed substantially due to the freedom of flipping from the 
^5*0 configuration to the ^Dq configuration and back. For A in the range of 1350 MeV-1600 MeV, we can find a bound 
state with the binding energy s = 6.46 — 168.73 MeV. The D wave probability increases with the binding energy e, 
and it is about 12.73% for e = 97.73 MeV, the importance of the tensor force is obvious. If all the coupling constants 
are reduced by 20%, we need increase A by about 200 MeV in order to obtain similar binding energy. However, the 
value of A is still in the reasonable range. Since the molecular state is widely extended, the decay into light mesons 
via annihilation is generally suppressed by the form factor. The leading source of decay is dissociation, to a good 
approximation the dissociation will proceed via the free space decay of the constituent mesons. Consequently the O"*""*" 
D*D* molecule mainly decays into DfD~jj, and Df DJjtt'^ , and the mode DJtt^tt'^ is forbidden by the phase 
space. 

For the axial vector 1"' state, there are two channels ^Si and ^Di. The coupling between the S wave and D wave 
has the same strength as the O"*""*" state, while the S wave spin-spin interaction potential Vs{r) is weaker than the 
corresponding one of the 0^^ state. Therefore the predictions for the binding energy and the static properties have 
similar pattern with the 0^+ sector, and the binding energy of the 1"' state is somewhat smaller than that of the 
latter for the same A value. We note that the unnatural spin parity forbids its decay into DsDs, while the decay 
mode DsD*/D*Ds is allowed. 

The state involves only one channel '^Pq. In contrast with the 0++ and 1^ cases, the tensor interaction 
potential — 2Vt(?') is attractive as a first order effect instead of a second order effect. The contributions of both 
spin-spin interaction and spin-orbit interaction are attractive as well, since Vs{r) and VLs{r) are positive as shown 
in Fig. [1] The potential Eq.© for this pseudoscalar is displayed in Fig. [5^ with A = 1600 MeV, we see that the 
potential is strong enough so that the P wave centrifugal barrier can be partly compensated, then there remains a 
weak attractive interaction in the intermediate range. Therefore bound state solution can be found for reasonable 
values of A, as can be seen from Table HVl For A = 1500 — 1600 MeV, we find the binding energy e = 1.40 — 114.81 
MeV. The binding energy is more sensitive to A than the 0++ and 1^ two coupled channels cases. 
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FIG. 2: The potentials for the single channel states with A = 1600 MeV. (a), (b) and (c) correspond to 0^"*", 1^^ and 1""*" 
states respectively. The solid line represents the potential from one boson exchange, and the dashed line denotes the effective 
potential after taking into account the centrifugal barrier. 



The results for 1++ state is similar to the single channel ^ case. Because the D wave centrifugal barrier is higher 
than the P wave centrifugal barrier, the total effective potential is less attractive than the ^ state, this point can 
be seen clearly in Fig. ^ If the coupling constants are reduced by 20%, bound state solution can be found only for A 
larger than 1920 MeV. 1"*"+ D*D* is harder to be bound than the ^ state due to the repulsive D wave centrifugal 
barrier. 

We then come to the CP exotic 1 and 2"' states, only one channel is involved in both sectors. As has been shown 
in Eq.® and Eq. pU)) . the potentials from one boson exchange are exactly the same, and they are less attractive than 
the potentials of the and 1++ states. For the 1 ^ configuration, bound state solution appears only for A as large 
as 2000 MeV, and we can find a 2"' bound state only if the regularization parameter A is larger than 3160 MeV. If 
we reduce all the coupling constants by 20%, larger value of A is required to find bound state solutions. Because the 
value of A is so large that it is far beyond the range of 0.8 to 1.5 GeV favored by the nuclcon-nucleon interactions, 
we tend to conclude that the CP exotic 1 ^ and 2"' D^Dl states can not be bound by the one boson exchange 
potential. This conclusion is consistent with the fact that no such CP exotic states have been observed so far. 

For the 1 states, there are three configurations ^Pi, ^Pi and ^Fi. In spite of the P wave centrifugal barrier, bound 
state solutions can be found for reasonable value of A, the reason is the large attractive contributions from the tensor 
interaction and the spin-orbit interaction. From the numerical results in Table IVIIIi we see that the binding energy 
is rather sensitive to A, and ^Pi is the dominant component. This is because that the 22 element of the potential in 
Eq. (fTTj) is more attractive than the 11 element, and the ^Fi component is strongly suppressed by the large F wave 
centrifugal barrier. The 1^^ D^D* state is very interesting, it can be produced via the e+e~ annihilation or with 
the help of the initial state radiation (ISR) technique at B factory. The existence of such a state can be confirmed 
or rejected, if more detailed e+e^ annihilation data in the range of 4100-4200 MeV become available. We strongly 
urge the Babar and Belle collaboration to search for this state, so that our prediction can be checked. In addition 
to the dominant decay modes and D'^D^jtt'-\ the 1 D*sD* molecule can also decay into DgDs and 

DsD*/D*Ds, which are other important decay modes. 

For the 2++ states, there are four channels ^D2; ^82, and ^G2- We note that both the tensor interaction and 
the spin-orbit interaction vanish in the ^82 configuration, and the spin-spin interaction is repulsive. However, bound 
state solution can be found for appropriate value of A albeit its value is somewhat larger than the corresponding ones 
in the 0+"'', 1"' , ^ and 1^~ cases. This is because that the mixing of with ^^2, ^£'2 and ^G2 under the 
tensor force increases the binding of the system through higher order iterative processes. From the numerical results 
in Table HXl we see that the ^^2 component is dominant, ^D2 probability is larger than ^D2, and ^G2 has the smallest 
probability. 

Finally two states 2 ^ and 2 remain. For both states, one has a two coupled channels (P wave and F wave 
configurations). Because of the P wave and F wave centrifugal barrier, bound state appears only for A as large as 
2080 MeV and 1890 MeV respectively. Therefore both the 2 ^ and 2 -D*Z)* states may not be bound by the one 
boson exchange potential. 

In short summary, ten allowed D^Dl states with low spin parity have been studied. We find that J^*^ = 0++, 1"' , 
'", 2++ and 1 states are more tightly bound. We expect that the 0+"*" and 1"*"" D^D* molecular states very likely 
exist, where the system can be in S wave or D wave and the two configurations mix with each other under the tensor 
interaction. Due to the remarkably strong contributions from the spin-spin interaction, the tensor interaction and the 
spin-orbit interaction, the D*Dl molecule may exist in spite of P wave centrifugal barrier. The 1~~ Dt^t state 
likely exist as well, and the ^Pi component is dominant. One expects that the one boson exchange potential could 
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support the 2++ D*D* bound state because of the presence of the S wave configuration and mixing with the other 
three liigher partial waves. However, the CP exotic 1 and 2"' states, 2 ^ and 2~~ states should not exist. 

IV. THE INTERPRETATION OF y(4140) AS A D'.D^ MOLECULE AND ITS IMPLICATION 

Since the one boson exchange potential does not support 1 '", 2^ , 2 and 2^^ D*Dl molecular states, we 
concentrate on 0++, 1^ , '", l"*"^, 1~~ and 2++ states in the following. Because the C— parity of y(4140) is 
positive, J^^ — 1^ and 1^^ states are not possible. From the numerical results in Table HB Tabic IXll it is obvious 
that the 0+"'' D*D* is most tightly bound by the one boson exchange potential. Consequently the most favorable 
quantum number of y(4140) is J^*^ = 0++, however, J^'~^ = or 2++ can not be excluded at present. It is crucial 
to perform a partial wave analysis in future, if the spatial parity turns out to be positive, J^^ = and 2++ are 
favored, otherwise it may be a state. 

If y(4140) is confirmed to be a 0"*"^ D^D* molecule with masses about 4143 MeV by future experiments, ^ and 
2++ D^D* molecules should be observed as well, whose masses are in the range of 4190 MeV to 4210 MeV. For the 
parameters that allow F(4140) to emerge as a 0+"*" D*D* molecule, one expects that the 1^ and 1 D*D* states 
may exist as well. We note that the C— parity of both states is negative. The 1~~ state is particularly interesting, 
its mass is about 4120-4150 MeV, and it can be produced directly in the e+e~ annihilation or via the initial state 
radiation at B factory. Both states can decay into DsD* / D*Ds, whereas the spin parity of the first state forbids its 
decay into DgDs- Since y(4140) is observed in the J/ip<f> channel, the 1^ and 1^^ states should be searched for in 
the J/tpii and J/ipr]' final states. We suggest the CDF, Babar and Belle collaboration to search for this state, which 
would be a critical check to the D*Dl molecule interpretation of y(4140). 

If F(4140) is identified to be a D*D* molecule, it is interesting to investigate whether a D*D* molecular state 
exists and what its most favorable quantum numbers are. In Ref. Q, the authors identified y(3940) as the D*D* 
molecular partner of F(4140), then its binding energy is about 77.5 MeV (Belle Collaboration) or 105.9 MeV (Babar 
Collaboration). However, in one pion exchange model Tornqvist demonstrated that the D*D* molecule should be 
near the threshold of 4020 MeV with J^<^ = 0++, 0"+, l+~ or 2++ It is necessary to reanalysis the D*D* system 
dynamically. In Ref. [ssj . the contributions from tt, t], cr, p and ui exchanges are included, the binding energy and other 
static properties are found to be sensitive to the regularization parameter A. However, we qualitatively confirmed 
that the J^*^ = 2++, 1 , 0++ and ^ D*D* states are more deeply bound. Consequently its huge binding energy 
implies the existence of more D*D* molecules with different quantum numbers, if F(3940) is identified as a D*D* 
bound state. However, no such candidates have been reported experimentally so far, therefore we tend to conclude that 
y(3940) as a D*D* molecule is not favored. Further experimental and theoretical efforts are needed to understand 
the structure and the properties of y(3940). With the same argument as that in the introduction section, F(3940) 
as a canonical cc(2P) charmonium can not be completely excluded. The charmonium assignment of y(3940) can be 
tested by searching for DD and DD* / D* D final states and by studying their angular distributions [3^ . 

For the B*B* system, the repulsive kinetic energy is greatly reduced due to larger mass of B* meson, therefore the 
B*B* state should be more deeply bound than D*D*. We expect that at least O"'"*' B*B* molecular state should exist 
with larger binding energy. Obviously this state could be searched for in the T{lS)(j) channel. Because of its large 
mass, the most promising places to produce this state conspicuously are the large hadron colliders such as Tevatron 
and LHC. 



V. CONCLUSIONS AND DISCUSSIONS 

In this work, we have dynamically studied the possible D^D* molecular states and the interpretation of F(4140) 
as a D*D* molecule in the one boson exchange model, where cr, rj and (j) exchanges are taken into account. Ten 
allowed states with low spin parity have been considered, we would like to stress that only S wave configuration is 
usually considered in the familiar phenomenological models such as the one boson exchange model in heavy quark 
effective theory [g^, HJ and the potential model with pairwise interactions [l^ [l^ . We find that the binding energy 
and static properties are sensitive to the regularization parameter A and the effective coupling constants. Since the 
regularization parameter A is poorly known so far, we are not able to precisely predict the binding energies of the 
possible D*Dl molecular states bound by one boson exchange potential. However, we can reliably predict which ones 
of the ten allowed states are much easier to be bound, and the predictions are rather stable even if the uncertainty 
of the coupling constants is considered, as is obvious from the numerical results listed in the manuscript. Further 
research on X(3872), which is a promising DD* / DD* molecule, would severely constraint the parameters of the one 
boson exchange model, especially the regularization parameter A, so that the predictions presented in the work could 
become more precise. 
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We quantitatively confirm that the 0+"'" D^D* state is most easily to be bound. Our numerical results imply that 
the J^^ = 0++, 0"+, 2++ and 1"" configurations are rather more strongly bound so that the corresponding 

D*D* molecules may exist, whereas the CP exotic 1 and 2^ , 2 ^ and 2^^ D*D* states are not be bound by 
the one boson exchange potential. We note that the possible existence of a number of bound state is not a specific 
prediction of our model (2lj, (23. Finally we would like to stress that we still can not completely rule out y(4140) as 
a conventional cc charmonium at present, in spite of its peculiar decay mode J/ilxj). From the theoretical predictions 
for the charmonium spectrum, F(4140) is most likely to be the 2^Z?2 state with I'^{J^'^') = 0+(2 '"), if it is a 
cc charmonium state. It unusual large branch ratio into J/tp4> may be explained by the rescattcring mechanism 
or the mixing between charmonium and molecule. Compared with other charmonium like states, the experimental 
information for y(4140) is scarce, and further experiment data are critically needed. 

For y(4140) as a D*D* molecule, we suggest that its most favorable quantum numbers are J^'~^ = although 
J-^*-^ = ^ and 2++ can not be ruled out by the present experimental data. It mainly decays into DfD~jj 
and D^DJjtt'^ via almost free decay of D* and D*, and the decay mode D^Dj ti^tt^ is forbidden by phase space 
constraints. The search for the four body decays y(4140) D'^Dj^^ and F(4140) DfD^jrr'^ is crucial to 
test the hadronic molecule hypothesis of y(4140). If y(4140) is confirmed to be a D*D* molecule by future 
theoretical and experimental efforts, the ^ and 2++ partners should exist with mass in the range of 4190-4210 MeV. 
We argue that the 1"' and 1 D^D* states with negative C— parity should be observed as well. The 1^^ states can 
be produced largely in the e^e^ annihilation or with the help the initial state radiation at B factory, and detailed 
e+e" annihilation data near 4100~4200 MeV are important to confirm or refute the existence of such state. Both 
the 1+" and D*Dl states can be searched for in the J/iprj and J/tprj' final states. We strongly urge the CDF, 
Babar and Belle Collaborations to search for these two negative C— parity states, which would be another important 
test to the molecular hypothesis of y(4140) and the reliability of our one boson exchange model. 

If we identify y(3940) as a D*D* molecule, its large binding energy requires the existence of more D*D* bound 
states, which has not been observed so far. Therefore the interpretation of y(3940) as a D*D* molecule may not 
be favored in our opinion. With the present experiment data, the charmonium assignment for F(3940) can not be 
ruled out. We suggest that the 0++ B*B* molecular state should exist, and it is bound more tightly than D*D*. We 
should search for this state at Tevatron or LHC in the T(1S')(/) channel. 
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APPENDIX A: NUMERICAL RESULTS FOR THE TEN ALLOWED DtD* STATES 



A(MeV) 


e(MeV) 


r,.ms(fni) 


Ps : Pd(%) 


1350 


6.46 


1.49 


99.48:0.52 


1400 


16.09 


1.03 


98.83:1.17 


1450 


32.11 


0.79 


97.51:2.49 


1500 


57.33 


0.63 


94.56:5.44 


1550 


97.73 


0.52 


87.27:12.73 


1600 


168.73 


0.44 


69.44:30.56 


all couplings are 


reduced by 20 percents 


A(MeV) 


e(MeV) 


rrms(fm) 


Ps : Pd(%) 


1500 


4.20 


1.75 


99.51:0.49 


1550 


10.28 


1.20 


98.94:1.06 


1600 


20.14 


0.91 


97.88:2.13 


1650 


35.21 


0.73 


95.77:4.22 


1700 


58.25 


0.60 


91.35:8.65 


1750 


95.26 


0.50 


81.68:18.32 



TABLE IL The predictions for the static properties of the J^'^' = 0"*"^ D*D* hadronic molecule, where e denotes the binding 
energy, rms is the root of mean square radius, Ps and Pd represent the S state and D state probabilities respectively. 



A(MeV) 


e(MeV) 


rrms(fm) 


Ps : Pd(%) 


1350 


6.13 


1.53 


99.65:0.35 


1400 


13.34 


1.12 


99.36:0.64 


1450 


23.76 


0.89 


98.95:1.05 


1500 


37.74 


0.75 


98.35:1.65 


1550 


55.74 


0.65 


97.48:2.52 


1600 


78.42 


0.58 


96.16:3.84 


1650 


106.87 


0.52 


94.15:5.85 


all couplings are 


reduced by 20 percents 


A(MeV) 


e(MeV) 


rrms(fm) 


Ps : Pd(%) 


1600 


10.55 


1.20 


99.29:0.71 


1650 


17.09 


0.98 


98.91:1.09 


1700 


25.48 


0.84 


98.39:1.61 


1750 


35.97 


0.74 


97.67:2.33 


1800 


48.92 


0.66 


96.67:3.33 


1850 


64.84 


0.59 


95.27:4.73 


1900 


84.49 


0.54 


93.29:6.71 


1950 


109.01 


0.49 


90.48:9.52 



TABLE III: The predictions about the binding energy, the root of mean square radius(rms) and the probabilities of the different 
components for the 1^^ D*D* molecule. 
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A(MeV) 


e(MeV) 


rrms (fm) 


1500 


1.40 


1.52 


1520 


14.99 


0.86 


1540 


33.36 


0.70 


1560 


56.13 


0.62 


1580 


83.26 


0.56 


1600 


114.81 


0.51 


all couplings are reduced by 20 percents 


A(MeV) 


M(MeV) 


rrms(fm) 


1700 


6.56 


0.99 


1720 


19.69 


0.76 


1740 


35.86 


0.65 


1760 


54.89 


0.58 


1780 


76.74 


0.54 


1800 


101.41 


0.50 



TABLE IV: The predictions for the binding energy and the rms of the ^ Dg Dg molecule. 



A(MeV) 


£(MeV) 


rrms(fm) 


1680 


4.21 


0.62 


1690 


28.16 


0.54 


1700 


54.83 


0.50 


1710 


84.04 


0.47 


1720 


115.71 


0.45 


all couplings are reduced by 20 percents 


A(MeV) 


e(MeV) 


rrms (fm) 


1920 


12.37 


0.54 


1930 


33.19 


0.49 


1940 


55.70 


0.46 


1950 


79.81 


0.44 


1960 


105.50 


0.43 



TABLE V: The predictions for the binding energy and the rms of the 1++ D*D| molecule. 



A(MeV) 


e(MeV) 


rrms(fm) 


2000 


10.72 


0.77 


2020 


23.60 


0.63 


2040 


38.60 


0.56 


2060 


55.61 


0.51 


2080 


74.60 


0.47 


2100 


95.54 


0.44 


2120 


118.44 


0.42 


all couplings are reduced by 20 percents 


A(MeV) 


e(MeV) 


rrms(fm) 


2480 


9.94 


0.72 


2500 


20.67 


0.60 


2520 


32.81 


0.53 


2540 


46.25 


0.49 


2560 


60.96 


0.45 


2580 


76.92 


0.42 


2600 


94.11 


0.40 


2620 


112.54 


0.38 



TABLE VI: The predictions for the binding energy and the rms of the 1 ^ DgDg molecule. 



A(MeV) 


e(MeV) 


rrms(fni) 


3160 


6.48 


0.37 


3170 


22.41 


0.34 


3180 


39.00 


0.33 


3190 


56.20 


0.32 


3200 


74.00 


0.31 


3210 


92.40 


0.30 


3220 


111.38 


0.29 


all couplings are reduced by half 


A(MeV) 


e(MeV) 


rrms (fm) 


4420 


11.51 


0.28 


4430 


28.97 


0.26 


4440 


46.92 


0.25 


4450 


65.36 


0.25 


4460 


84.27 


0.24 


4470 


103.63 


0.23 



TABLE VIL The predictions for the binding energy and the rms of the 2^ DgDg molecule. 
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A(MeV) 


e(MeV) 


rrms(fm) 


Ppo : Pp2 : Pp(%) 


1470 


5.38 


1.10 


3.40:96.50:0.10 


1480 


13.13 


0.89 


3.18:96.71:0.11 


1490 


22.23 


0.78 


3.01:96.87:0.13 


1500 


32.59 


0.71 


2.86:97.00:0.14 


1510 


44.18 


0.66 


2.73:97.12:0.16 


1520 


56.97 


0.62 


2.61:97.22:0.17 


1530 


70.95 


0.59 


2.51:97.30:0.19 


1540 


86.12 


0.56 


2.42:97.38:0.20 


1550 


102.49 


0.53 


2.34:97.44:0.22 


all couplings are 


reduced by 20 percents 


A(MeV) 


£(MeV) 


rrms(fm) 


Pp():Pp2:Pf(%) 


1650 


11.28 


0.88 


2.46:97.45:0.09 


1660 


18.49 


0.78 


2.36:97.53:0.11 


1670 


26.53 


0.71 


2.29:97.60:0.12 


1680 


35.36 


0.67 


2.22:97.66:0.13 


1690 


44.96 


0.63 


2.15:97.71:0.14 


1700 


55.33 


0.60 


2.10:97.75:0.15 


1710 


66.45 


0.57 


2.04:97.79:0.16 


1720 


78.32 


0.55 


2.00:97.83:0.18 


1730 


90.95 


0.53 


1.95:97.86:0.19 


1740 


104.33 


0.51 


1.91:97.89:0.20 



TABLE VIII: The predictions about the binding energy, the root of mean square radius and the probabilities of the different 
components for the 1~~ D*D* molecule, where Ppo and Pp2 denote the ^Pi state and ^Pi state probabilities respectively. 



A(MeV) 


£(MeV) 


ri.nis(fm) 


Pdo :Ps :Pd2 :Pg(%) 


1350 


6.81 


1.47 


0.07:99.48:0.46:0.00 


1400 


12.03 


1.18 


0.10:99.14:0.75:0.00 


1450 


18.59 


1.00 


0.15:98.67:1.18:0.00 


1500 


26.55 


0.88 


0.21:98.01:1.78:0.00 


1550 


36.10 


0.80 


0.29:97.07:2.63:0.00 


1600 


47.59 


0.73 


0.40:95.74:3.85:0.01 


1650 


61.65 


0.67 


0.53:93.85:5.60:0.02 


1700 


79.28 


0.62 


0.70:91.12:8.13:0.05 


1750 


102.10 


0.57 


0.91:87.18:11.77:0.14 


all coupling 


;s are reduced by 20 percents 


A(MeV) 


e(MeV) 


ri.nis(fm) 


Pdo :Ps :Pd2 :Pg(%) 


1650 


8.09 


1.36 


0.11:98.92:0.97:0.00 


1700 


11.69 


1.18 


0.15:98.44:1.40:0.00 


1750 


16.09 


1.05 


0.21:97.81:1.98:0.00 


1800 


21.41 


0.94 


0.27:96.97:2.75:0.01 


1850 


27.91 


0.86 


0.35:95.86:3.79:0.01 


1900 


35.91 


0.78 


0.44:94.36:5.17:0.02 


1950 


45.91 


0.72 


0.56:92.36:7.04:0.04 


2000 


58.62 


0.66 


0.70:89.64:9.56:0.09 


2050 


75.18 


0.60 


0.89:85.95:12.95:0.22 


2100 


97.42 


0.55 


1.13:80.82:17.48:0.57 


2150 


129.26 


0.50 


1.49:73.08:23.51:1.93 



TABLE IX: The predictions about the binding energy, the root of mean square radius and the probabilities of the different 
components for the 2^+ DgDs molecule. 
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A(MeV) 


e(MeV) 


rrms(fm) 


Pp : Pf(%) 


2080 


21.35 


0.54 


32.02:67.98 


2090 


50.76 


0.46 


25.30:74.70 


2100 


84.10 


0.42 


21.08:78.92 


2110 


120.80 


0.39 


18.09:81.91 


all couplings are 


reduced by 20 percents 


A(MeV) 


e(MeV) 


r,.ms(fm) 


Pp : Pf(%) 


2510 


18.99 


0.47 


24.99:75.01 


2520 


47.32 


0.40 


20.26:79.74 


2530 


78.24 


0.37 


17.38:82.62 


2540 


111.36 


0.35 


15.34:84.66 



TABLE X: The predictions about the binding energy, the root of mean square radius and the probabilities of the different 
components for the 2""*" DgDg molecule. 



A(MeV) 


e(MeV) 


ri.ms(fm) 


Pp :Pf(%) 


1890 


1.59 


1.12 


65.79:34.21 


1900 


21.17 


0.61 


51.31:48.69 


1910 


46.97 


0.51 


41.55:58.46 


1920 


78.03 


0.46 


34.08:65.92 


1930 


113.83 


0.42 


28.31:71.69 


all couplings are 


reduced by 20 percents 


A(MeV) 


e(MeV) 


r,.ms(fm) 


Pp :Pf(%) 


2230 


5.92 


0.71 


49.37:50.63 


2240 


28.13 


0.50 


38.73:61.27 


2250 


54.53 


0.44 


31.95:68.05 


2260 


84.34 


0.40 


26.99:73.01 


2270 


117.14 


0.38 


23.20:76.80 



TABLE XI: The predictions about the binding energy, the root of mean square radius and the probabilities of the different 
components for the 2~~ D*D* molecule. 



